Fingerprinting Dark Energy II: weak lensing and galaxy clustering tests. 
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The characterization of dark energy is a central task of cosmology. To go beyond a cosmological 
constant, we need to introduce at least an equation of state and a sound speed and consider observa- 
tional tests that involve perturbations. If dark energy is not completely homogeneous on observable 
scales then the Poisson equation is modified and dark matter clustering is directly affected. One 
can then search for observational effects of dark energy clustering using dark matter as a probe. In 
this paper we exploit an analytical approximate solution of the perturbation equations in a general 
dark energy cosmology to analyze the performance of next-decade large scale surveys in constrain- 
ing equation of state and sound speed. We find that tomographic weak lensing and galaxy redshift 
surveys can constrain the sound speed of the dark energy only if the latter is small, of the order of 
c s < 0.01 (in units of c). For larger sound speeds the error grows to 100% and more. We conclude 
that large scale structure observations contain very little information about the perturbations in 
canonical scalar field models with a sound speed of unity. Nevertheless, they are able to detect the 
presence of "cold" dark energy, i.e. a dark energy with non-relativistic speed of sound. 

PACS numbers: 98.80.-k; 95.36.-fx 
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I. INTRODUCTION 

Even though there is ample observational evidence that 
the expansion rate of the Universe is accelerating [l|, [2| , 
there is still no convincing theoretical model that can ex- 
plain the observations. The simplest model in agreement 
with the data is the cosmological constant. However, it 
suffers from size and coincidence problems. In addition, 
it is likely that an earlier phase of accelerated expansion 
happened that was not due to a cosmological constant 3] . 
This period, inflation, is normally modeled as being due 
to a scalar field. The same basic mechanism can explain 
the dark_energy, and is often called Quintessence in this 
context [4]. Scalar field dark energy is often employed 
to model a general equation of state w{z), although one 
can also build a perfect fluid with exactly the same first 
order perturbations [5] . 

Over the last few years it has become clear that w alone 
is not sufficient to distinguish between different models 
explaining the dark energy, and that the perturbations 
induced by the dark energy constitute a complementary 
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probe |fML2|. The contribution of dark energy fluctu- 
ations to the perturbation dynamics is generally small 
and can be very difficult to measure. In a recent paper 
some of us derived theoretical approximations to the per- 
turbations generated in scalar field dark energy models 
13]. In this paper, we will use these expressions to inves- 
tigate how two cosmological probes that will be of great 
importance in the future, weak lensing (WL) and the full 
galaxy power spectrum, can be used to measure the typ- 
ical features that appear in such models, specifically the 
existence of a sound horizon. Exploiting the advantage of 
analytical approximations, we will try to isolate in some 
detail where the signal comes from. 

Anticipating the conclusion, we will show that even 
a full sky, tomographic, deep redshift and imaging sur- 
vey cannot constrain significantly the sound speed of the 
dark energy unless it is smaller than, roughly, 0.01c 1 . 
Small values of c s are not compatible with a scalar field 
description with standard kinetic term and require ex- 
tensions like k-essence 14j but are in general not forbid- 
den: we will show that even with c s — 0.01, dark energy, 
although "cold" , remains much less clustered than dark 
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1 Here c denotes the speed of light, and in the rest of the paper we 
will use units so that c = 1. 



matter and will not develop non-linearities by the present 
time. The special case of c s = has some additional at- 
tractive features, in that such models can cross w = — 1 
p, [l5[ and can act as unified models of dark energy and 
dark matte r 1161) (but are also subject to the "dark de- 
generacy" [17|), For such extreme cases however dark 
energy becomes non linear along with dark matter and 
appropriate non-linear corrections are needed, especially 
for the weak lensing probe. 

The sound speed of the dark energy could be mea- 
sured also by other effects, like the Integrated Sachs- 
Wolfe (ISW) effect, either on the CMB or cross-correlated 
with large-scale structure. In [18| it has been shown that 
the CMB ISW tail can put a lower limit to c 2 , while [l9| 
finds that c s « 0.01 is actually favored by the WMAP 
3-year data. The ISW is complementary to the probes 
we consider here since it is sensitive to perturbations on 
larger scales. 

As in [13[ we set the anisotropic stress of the dark en- 
ergy to zero, which is appropriate for scalar field dark en- 
ergy and in general for standard gravity models. In order 
to investigate more general cases, for example modifica- 
tions of General Relativity, a non-zero anisotropic stress 
would have to be taken into account [10|, but this is left 
for future work. 

In detail, the paper is organized as follows. We be- 
gin by discussing the main results obtained in [1J|. We 
then use these results to study the impact of the scalar 
field dark energy perturbations on the weak lensing sig- 
nal and on probes of galaxy clustering, paying special 
attention to whether the presence of a sound horizon can 
be detected (by measuring the sound speed c 2 ). We also 
discuss which aspect of these probes are most sensitive 
to the perturbations. 

We finally investigate constraints on the integrated de- 
viation from the Poisson equation, which we call W, 
which may be useful to quantify the power of probes be- 
yond the conventional, w-based Figure of Merit used by 
the dark energy task force j20j | . 



II. SETTING THE SCENE 

In this section we define our notation and present a 
short discussion of the perturbation equations and the 
approximate analytical solution of [13]. We refer the 
reader to that paper (and citations therein) for more de- 
tails. 



Hubble parameter that we will consider here is: 
H 2 = (^\ = Hi [f! m ,oa- 3 + (1 - On,o) 9(a)] (1) 



do 



and the subscript 



where g(a) — exp —3 J — ^£W. 

denotes the present epoch. This expression implies that 
we limit ourselves to a fiat universe which is filled with 
matter (designated by a subscript m) and a general fluid 
with average pressure p = wp. In general the equation of 
state parameter w is a function of time, but we will take 
it to be constant later on. 

We will consider linear perturbations about this 
spatially-flat background model, defined by the line el- 
ement: 

ds 2 = a 2 [- (1 + 2V0 dr 2 + (1 - 20) d^dz 1 ] . (2) 

As is apparent from the line element, we use the con- 
formal Newtonian (longitudinal) gauge and retain only 
scalar perturbations. The perturbation equations for a 
dark energy fluid with sound speed c s and parameter of 
state w are: 
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where V — (1 + w)(ikiV l ) is a measure of the velocity 
perturbation, 5 is the density contrast, and the prime 
here means the derivative with respect to the scale factor 
a. We parameterize the pressure perturbation as: 
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where c 2 



w 



is the adiabatic sound speed 



3H{l+w) 

and c 2 is the sound speed in the rest-frame of the dark 
energy fluid. As already mentioned, we assume that the 
parameter of equation of state w stays constant so that 
c 2 =w. 

We also assume a vanishing anisotropic stress a = 
as is the case for Quintessence and K-essence models, 
and therefore we have ip = <f>. Finally, the gravitational 
potential can be found with the help of the Einstein equa- 
tions, 
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A. Definitions 

In the following, the dots will refer to the derivatives 
with respect to the conformal time r which is related 
to the universal time t by dt = a (t) dr. The physical 



B. Analytical solutions for dark energy 
perturbations 

If the universe is perfectly matter dominated, then k 2 <p 
is a constant. In [13| we used this result to derive ana- 
lytical solutions to the perturbation equations for dark 
energy in different limits: 



• scales larger than the sound horizon. The 

modes larger than the sound horizon are found for 
k <C aH/c s . In this case, we neglect all terms con- 
taining the sound speed in Eq. @ and then in 
Eq. <j3j> , effectively setting c 2 = 0. Neglecting a 
decaying solution for which V ex l/a we find: 
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Scales smaller than the sound horizon. For 

modes smaller than the sound horizon, k ^> aH/c s , 
we consider only the terms containing k 2 . We 
find that the density perturbations become con- 
stant and that the velocity perturbations decay: 
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In these equations, the factor 8i n sets the overall scale of 
the perturbations, chosen so that the matter perturba- 
tions are 8 m — 8i n a on sub-horizon scales. In addition, 
the dark energy perturbations depend on their present 
density fraction (equal to 1 — Q, m ^ since the universe 
is taken to be flat) but only through the combination 
Qm,oh 2 , the equation of state parameter w and the sound 
speed c 2 . They are a function of both scale k and scale 
factor a. 

Let us quickly compare dark energy and dark matter 
fluctuations. For large scales the dominant term in Eq. 
(j7)l is the term containing k. Dark matter has a similar 
solution except that wdm = 0; the ratio of the squares 
(i.e. the power spectrum ratio) for large scales then will 
be : 
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if w < —0.8. However, for small scales (but still larger 
than c s /aH), the first term in Eq. (|7|) dominates; in this 
case the ratio is: 
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0.0035. 



(12) 



The dark energy perturbations are thus always sup- 
pressed relative to the dark matter perturbations, even 
for a very low sound speed (note that the solution Eq. 
([7|) was obtained assuming c 2 =0 ). 

For the forecasts of this paper, it is important to know 
whether we have to worry about non-linear effects: the 
non-linear clustering of dark energy models is still badly 
understood and may depend on the precise action of the 
model, and it is not clear what corrections we would 



need to apply. Outside the sound horizon, the dark en- 
ergy perturbation track those in the dark matter by the 
factors given above. These are roughly upper limits as 
w = —0.8 is at the upper limit of what is still allowed, 
and values of w closer to — 1 will lead to smaller dark en- 
ergy perturbations. If c s = 0.01 the sound horizon at the 
present is of the order of c s (2/ Ho) f=a 60 Mpc/h. On these 
scales dark matter fluctuations are still linear and there- 
fore dark energy fluctuations are safely below the non- 
linearity threshold. On sub-sound-horizon scales, dark 
energy perturbations stop growing and become quickly 
negligible relative to the dark matter perturbations. If 
c s is smaller than 0.001 then the dark energy will behave 
much like dark matter on the cluster scales and an appro- 
priate non-linear correction would be needed, especially 
on the small scales probed by weak lensing (the small 
scale cut-off in the galaxy clustering safely excludes non 
linear scales). Therefore we do not extend our analy- 
sis to such small sound speeds; we include however for 
completeness the extreme case c 8 = 0, for which we as- 
sume that we can apply the same non-linear correction 
as for the dark matter - whether this is acceptable may 
in general depend on the precise action of the field. 



III. IMPACT OF THE SOUND SPEED ON 
GALAXY CLUSTERING AND WEAK LENSING 

In this section we discuss the effect of dark energy clus- 
tering on the galaxy and weak lensing probes. We derive 
approximate analytical expressions that show how the 
relevant quantities change with the sound speed. This 
will help understanding qualitatively the numerical re- 
sults of the next section. 



A. The Q parameter 

Matter domination was a necessary ingredient to derive 
the solutions given above. However, dark energy comes 
to dominate eventually, and then the potential starts to 
decay and the perturbations grow more slowly or start 
to decrease. 

It is difficult to capture this behavior accurately. A 
way around this problem can be found by looking at the 
variable Q(k, a) which we introduced in [12[ to describe 
the change of the gravitational potential due to the dark 
energy perturbations. Q is defined through 
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If the dark energy or modification of gravity does not 
contribute to the gravitational potential (for example if 
the dark energy is a cosmological constant) then Q = 1. 
Otherwise Q will deviate from unity, and in general it is 
a function of both scale and time. 



Introducing the comoving density perturbation A 
5 + 3aHV/k 2 , we can write 
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We neglect radiation in this paper as we are interested in 
observational tests at late times where it is subdominant. 
Then Q is defined as 
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Using simply the solution for the perturbations during 
matter domination, as discussed in the previous para- 
graph, we find that the resulting expression for Q is sur- 
prisingly accurate even at late times. The reason is that 
both fluids, dark energy and matter, respond similarly 
to the change in the expansion rate so that most of the 
deviations cancel. We find that the sub-soundhorizon ex- 
pression below is accurate at the percent level, while on 
larger scales there are deviations of about 10 to 20% by 
today (depending on w). The latter can be corrected "by 
hand" in order to obtain a more precise formula, but the 
expressions are sufficiently accurate for our purposes and 
we keep them as they are. 

As we have set the anisotropic stress to zero, the per- 
turbations are fully described by Q. In [13| we provided 
the following explicit expression for the Q (k, a) which 
captures the behavior for both limits (above and below 
the sound horizon): 



Q{k,a) = 1 + 
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Here we used v{a) 2 = k 2 c 2 s a/ (fi m ,o-ffo) which we defined 
through c s k = vaH so that v counts how deep a mode 
is inside the sound horizon. We see that during matter 
domination (used in the first expression for v) a k mode 
moves ever more deeply inside the sound horizon as a 
grows, which is of course no surprise. 

From Eq. (fT5|) we see that the deviation of Q from 1 
depends on the ratio of the energy density of the dark 
matter and the dark energy, and on the relative amount 
of perturbations. The former scales as a~ 3w for a con- 
stant w and the latter behaves as discussed in the last 
section. For fi m o = 0.25 and w = —0.8 we have that 
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on scales that are larger than the sound horizon, v ss 0. 
This is not a negligible deviation today, but it decreases 
rapidly as we move into the past, as the dark energy 
becomes less important. 2 As a scale enters the sound 
horizon, Q — 1 grows with one power of the scale factor 



2 For this reason, early dark energy models can have a much 
stronger impact |21| . 



slower (since Sde stops growing), suppressing the final 
deviation roughly by the ratio of horizon size to the scale 
of interest. In the observable range, (k/Ho) 2 « 10 2 — 10 4 . 
Therefore if c s « 1, Q — > 1 and the dependence on c s is 
lost. This shows that Q is sensitive to c s only for small 
values, c 2 < 10~ 2 . 

We can characterize the dependence of Q on the main 
perturbation parameter c 2 by looking at its derivative, a 
key quantity for Fisher matrix forecasts: 



<91ogQ 
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0-1 



(1 + x) 
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where x = |i/(a) 2 /(l - 3w) ~ 0.2v(a) 2 (with the last 
expression being for w = —0.8). For the values we are 
interested here, this derivative has a peak just inside the 
sound horizon, the exact position is 
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Today the sound horizon is given by c s ss Ho/k. It lies 
in the observable range of k for sound speeds of the order 
of c s s=s 0.01 — 0.001. We plot the derivative as the red 
curve in Fig. [51 there and from Eq. (fT8")) we can see that 
the shape is basically proportional to —x/(l + x) 2 . The 
derivative with respect to c 2 (instead of log c 2 ) contains 
an additional factor of 1/c 2 which means that it will be 
boosted by several orders of magnitude for small sound 
speeds. 

We will later forecast how well the deviation of Q from 
1 due to the dark energy perturbations can be measured 
by future cosmological surveys. 



B. The growth rate and the 7 parameter 

In the ACDM model of cosmology, the dark matter 
perturbations on sub-horizon scales grow linearly with 
the scale factor a during matter domination. During ra- 
diation domination they grow logarithmically, and also 
at late times, when dark energy starts to dominate, their 
growth is suppressed. It is well known that in ACDM the 
growth factor can be expressed as: 



G(a) 



_. S m (a) 
S m (ao) 



exp 



Sl m (a')' 



da 



'} (20) 



where 7 ~ 0.545 is called the growth index. There are 
two ways to influence the growth factor: firstly at back- 
ground level, with a different Hubble expansion. Sec- 
ondly at perturbation level: if dark energy clusters then 
the gravitational potential changes because of the Poisson 
equation, and this will also affect the growth rate of dark 
matter. All these effects can be included in the growth 
index 7 and we therefore expect that 7 is a function of 
w and c 2 (or equivalently of w and Q). 



According to [221 ], the growth index depends on dark 
energy perturbations (through Q) as 

3(l-w-A(Q)) 
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where 



A(Q) 
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However, dark energy perturbations (hence Q) are dif- 
ficult to measure at least if the dark energy has a very 
low sound speed. The growth index factor seems to be 
a more promising parameter and several experiments are 
planned to measure this quantity. Furthermore, we can 
invert Eq. (|2"Tj) and ask the question: in the absence of 
anisotropic stress which Q is needed to generate a given 
7? Using Eqs. dHJ and ^^ we have: 



Q - 1 = [1 - n m (a)} 



1 — w — — (5 — 6w) 7 



(23) 



The last equation is worth another look: let us assume 
we measure the growth index 7 = 6/11 (the value usually 
associated to the cosmological constant) then Eq. (|2"3")l 
becomes: 

Q - 1 = i (1 + w) [1 - n m (a)] . (24) 

If we are dealing with the cosmological constant then we 
have the expected result of Q— 1 = 0. However, if w 7^ —1 
then we can evaluate the corresponding value of Q; for 
instance, if w — —0.8 then we have Q — 1 ~ 1CP 2 , which 
is a positive number. However, this value is fairly big for 
dark energy perturbations, see [13| . Moreover, for non- 
phantom models the quantity Q — 1 is always positive 
due to the relative increase of dark energy perturbations 
and the (1 + w) factor; let us remind the reader that if 
the anisotropic stress is set to zero then Q fully describes 
the evolution of dark energy perturbations; so, we can 
think to rephrase the sentence as: if Q — 1 > which 7 
can we reach? Setting Eq. ([23]) positive we have 
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5 — 6w 

which sets a sort of upper bound to 7; moreover, this is 
the limit we have if dark energy perturbations are set to 
zero (which can be seen directly from Eq. (12TT0 . 

In general we can say that if the anisotropic stress is 
zero then dark energy perturbations always decrease the 
value of the growth index. If we want to allow for per- 
turbations in the dark energy sector while increasing at 
the same time the growth index then we need a non zero 
anisotropic contribution; this is the case for the DGP 
model (treated as an effective dark energy perturbations, 
see [l(|) where 7 ~ 0.68. 

Let us also consider the derivative of log G with respect 
c 2 to gain an idea of how strongly the growth factor de- 
pends on sound speed, and thus on the characteristics of 
the dark energy perturbations: 
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a dc 2 



log ft m (a) Q m (a) 7 da. (26) 



We can see more clearly the dependence on the param- 
eters by noting that the quantity log Q m (a) Q m (a) 7 can 
be approximated as J7 m (a) — 1 for Q m close to unity 
(we use however the exact expressions in the calculations 
below). Eq. (f26|) then reads: 
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where b — x/a ~ Q.2v 2 /a is a constant during matter 
domination (the value of x today). We first consider two 
different limits 

• 5>1 (sub-sound horizon regime): 

d log G 3 1 + w 1 — il m 
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• b ^C 1 (super-sound horizon regime) : 
d log G 3 1 + w 1 — tt m 
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Both of these derivatives, like the one of Q discussed 
further above, have a strong dependence on the sound 
speed. As before, the peak location is near the sound 
horizon, and again the derivative with respect to c 2 adds 
an additional 1/c 2 factor so that small sound speeds will 
lead to much larger values of the derivative. 

The numerical factor (including the w dependence) is 
—0.22 and —0.09 for w — —0.8, respectively. We also 
notice that the first limit does not depend on Q — 1 like 
the derivative Q, but on Q — Qq as it is an integral. The 
growth factor is thus not directly probing the deviation 
of Q from unity, but rather how Q evolves over time. 
This is similar in the second limit, except that we probe 
the evolution of (Q — l)x. 

We can find a unified formula which accounts for both 
regimes: 
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r 



where Qq and xo are evaluated at a = o,q = 1. 

The derivative is shown as the blue dashed line in 
Fig. [5] Just as the Q derivative, the G derivative peaks 
close to the sound horizon, but is smaller than the one 
of Q at this redshift. It grows however with growing red- 
shift due to its integral nature, while the other derivatives 
decrease. 
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C. Shape of the dark matter power spectrum 



In order to quantify the impact of the sound speed on 
the matter power spectrum we proceed in two ways: we 
first use the CAMB output [23J (which contains the full 
information on dark energy perturbations) and then we 
consider the analytic expression from Eisenstein & Hu 
[24( | (which does not include dark energy perturbations, 
i.e. does not include c s ). 

As anticipated, we find here that the impact of the 
derivative of the matter power spectrum with respect 
the sound speed on the final errors is only relevant if 
high values of c 2 are considered; for instance if c~ = 1 
then the errors on the sound speed increase by about 10 
times when the analytic formula in Ref. 2A\ is consid- 
ered instead of the CAMB output. By decreasing the 
sound speed, the results are less and less affected, and 
are completely unchanged when c 2 = 10~ 5 . The reason 
is that for low values of the sound speed other param- 
eters, like the growth factor, start to be the dominant 
source of information on c 2 , see Fig. [71 Since for large 
c s our results show that the sound speed in practically 
not measurable with the probes here considered, we dis- 
regard this difference between the analytical fit and the 
CAMB output and use the analytical fit whenever we 
find it convenient, in particular in estimating the weak 
lensing Fisher matrix. 

In Fig. [1] we show the derivative of the logarithm of 
the matter power spectrum (from CAMB) with respect 
to the logarithm of sound speed. It is easy to understand 
this result: In Fig. 3 of [13] we showed that the power 
spectrum decreases by a few percent in a smooth step- 
like fashion at a scale just inside the sound horizon. Here 
we see the derivative of this transition. Just like the step, 
the peak of the derivative is situated a little bit inside the 
sound horizon, as is the case for the other contributions 
as well. 



Figure 1: The figure shows the derivatives of the matter power 
spectrum with respect to log c 2 for two different values of 
the sound speed: c 2 — 1CP 2 and c 2 — 10~ 4 , red solid and 
blue dashed line, respectively. The vertical dotted lines give 
the scale of the sound horizon for two different sound speeds, 
c 2 = 1(T 4 (left) and c 2 = 1(T 2 (right). 



D. Redshift space distortions 

The velocity field of matter is directly sourced by 
the gradient of the tp gravitational potential, and as- 
tronomers have been trying to measure it for a long 
time [25|, [26[ . More recently constraints on the velocity 
perturbations from redshift space distortions became an 
important probe for dark energy phenomenology, since 
their combination with weak lensing allows to disentan- 
gle the two potentials, which in turn allows to put limits 
on the anisotropic stress of the dark sector [12| • This is 
an important test for modifications of gravity |1(|. Al- 
though there is no anisotropic stress in our scenario, we 
are nonetheless interested in how redshift space distor- 
tions are modified by the presence of dark energy pertur- 
bations. 

The distortion induced by redshift can be expressed in 
linear theory by the f3 factor, related to the bias factor 
and the growth rate via: 
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The derivative of the redshift distortion parameter with 
respect to the sound speed is: 

cnog(l + /V) 



d log c 2 
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We see that the behavior versus c 2 is similar to the one 
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Figure 2: The figure shows the derivatives with respect to 
logcj of the Q parameter Eq. (|18p (red solid line), the growth 
factor Eq. (I30|l (blue dashed line), the redshift distortion pa- 
rameter Eq. (|32p (green dotted line) and the matter power 
spectrum Po(fc) (black dot-dashed line). The derivatives are 
evaluated at a z = 0.5 except the one for Po(k) which is at 
2 = 0; for all the derivatives the sound speed is c? s — 10~ 4 . 
The vertical dotted line gives the scale of the sound horizon 



for ci 



10 -4 at z = 0.5. 



for the Q derivative, so the same discussion applies. Once 
again, the effect is maximized for small c s . 

In Fig. [5] we compare the derivative of the redshift 
space distortion (green dotted line, averaged over /x) with 
those of the growth factor G and the Q variable. The j3 
derivative is comparable to G at z = but becomes more 
important at low redshifts and, as mentioned, is similar 
to the Q derivative in shape. 



IV. OBSERVATIONAL CONSTRAINTS 
A. Weak lensing 

We want to investigate now the response of weak lens- 
ing (WL) to the dark energy parameters. We proceed 
with a Fisher matrix as in [12J ] , to which we refer for the 
implementation details, the main difference here being 
that the parameter Q, which was general, now has an 
explicit form that describes how it changes as a function 
of scale and redshift. Since Q depends on w and c 2 , we 
can forecast the precision with which those parameters 
can be extracted. We can also try to trace where the 
constraints come from. 

For WL experiments the important quantity is the 
lensing potential which is, using the metric Eq. ^: 



$ = <fr + ip. 



(33) 



For a vanishing anisotropic stress the WL potential be- 
comes: 



k 2 $ 



-2Q 



3H Q mi o 
2a 



A, ; 



(34) 



decomposed as: 

A m (a, k) = aG (a, k) A m (fc) . 



(35) 



Here A m (fc) = A m (a = l,fc) determines the matter 
power spectrum today, P(fc) = |A m (fc)| 2 and G(a,k) is 
the growth factor, which depends here not only time but 
also on fc since the dark energy perturbations in general 
introduce a scale dependence. 
We can write Eq. (f34|) as: 

fc 2 $ = -3H (a) 2 a 3 Q (a, k) Q, m (a) G (a, k) A m (fc) (36) 



In linear perturbation theory all fc modes evolve indepen- 
dently, so that the dark matter density contrast is often 



where we used Eq. (f35|). 

Hence, the lensing potential contains three conceptu- 
ally different contributions from the dark energy pertur- 
bations: 

• The direct contribution of the perturbations to the 
gravitational potential through the factor Q, see 
section MI Al 

• The impact of the dark energy perturbations on 
the growth rate of the dark matter perturbations, 
affecting the time dependence of A m , through 
G (a, fc), cf section ITlIB I 

• A change in the shape of the matter power spec- 
trum P(fc), corresponding to the dark energy in- 
duced fc dependence of A m , as discussed in section 
iHTCl 

We consider a typical next-generation tomographic weak 
lensing survey characterized by the sky fraction f s k y = 
1/2 and by the number of sources per arcmin 2 , d = 40. 
We consider the range 10 < I < 10000 and we extend 
our survey up to three different redshifts: z max = 2,3,4. 
For the non linear correction we use the halo model by 
Smith et al. [27| . We choose as fiducial model Q m o = 
0.24, h = 0.7, n DE = 0.737, Q K = 0, n b h 2 = 0.0223, 
r = 0.092, n s — 0.96, wq = —0.8, and several values for 
c 2 . The observationally borderline fiducial value of wq is 
chosen so as to maximize the impact on Q: values closer 
to —1 reduce the effect and therefore increase the errors 
on c s . 

We always plot the fully marginalized confidence el- 
lipses (and quote fully marginalized errors) at 68%, which 
in 2D correspond to semi-axes of length 1.51 time the 
eigenvalues. In Fig. [3] we report the confidence region 
for wq, c 2 for two different values of the sound speed and 
Zmax- For high value of the sound speed (c 2 = 1) we 
find cr(wo) = 0.0171 and the relative error for the sound 
speed is a(c 2 )/cl — 2768. As expected, WL is totally in- 
sensitive to the clustering properties of quintessence dark 
energy models when the sound speed is equal to 1. The 
presence of dark energy perturbations leaves a w and c 2 
dependent signature in the evolution of the gravitational 
potentials through A^£;/A m and, as already mentioned, 
the increase of the c 2 enhances the suppression of dark 
energy perturbations which brings Q — > 1. 
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Table I: Here are listed the errors for the equation of state 
parameter and the relative errors for the sound speed at 



Once we decrease the sound speed then dark energy 
perturbations are free to grow at smaller scales. As an ex- 
ample we show in the lower panel of Fig. [3] the confidence 
region for w$,c 2 s for c 2 s = 10~ 5 , we find a(wo) = 0.025, 
a(cj.)/c 2 = 0.49; in the last case the error on the mea- 
surement on the sound speed reduced to the 50% of the 
total signal. 

In Tab. U we list the errors for wq and the relative 
errors of c 2 for six different values of the sound speed for 
a survey up to z max = 3. 

We should note here that the Fisher matrix approach 
breaks down for such large errors since it only considers 
the local curvature of the likelihood at the location of the 
fiducial model 28]. If we performed an actual analysis 
for data with c 2 = 1, wc would find a lower limit near 
c 2 — 10 -5 since those models are measurably different 
from c 2 s — 1 according to the lower panel of Figure |31 
and the ellipse would be cut off there. On the other 
hand, there would not be any upper limit to c 2 s since 
for all those cases the dark energy perturbations are not 
detected. 

To explore more directly which term influences the 
most the weak lensing signal we assume that there are 
two different Q's: one, Q 7 , which enters directly in the 
growth index expression (|2ip and the other, Q$, which 
is the term that enters linearly in the expression of the 
gravitational potential (j34|) . We concentrate here only 
on these terms since, as we mentioned previously, the 
changes on the matter power spectrum are negligible in 
the interesting range (i.e. for small sound speed). 

• Q 7 only: Here we consider only the contribution 
from Q 7 setting effectively Q$ = 1. 

• Q$ only: Here we assume that dark energy pertur- 
bations enter only through Q$, effectively setting 
the function Q 1 = 1 in Eq. (|2"Tj) . 

We show the behavior of the errors in Fig. [4] as a func- 
tion of the sound speed c 2 s . We notice that the main 
contribution comes from Q$. Keeping only this term 
leads to the same errors as the full expression. This is 
consistent with the discussion in sections MI Al and IIIIBI 
and the derivatives shown in Fig. [5] where we see that the 
growth factor is always less sensitive to the dark energy 
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Figure 3: Confidence region at 68% for three different value 
of z max = 2, 3, 4, blue dashed, green long-dashed and solid 
contour, respectively. The upper panel shows the confidence 
region when the sound speed is c s = 1; the bottom panel with 
the sound speed <? s = 10~ 5 . The parameter equation of state 
is for both cases wq = —0.8. 




Figure 4: The figure shows the relative errors on the sound 
speed and their corresponding fits for our benchmark surveys; 
assuming only Q 7 in red and assuming only Q$ in blue. 



perturbations than Q. As both enter in the same way in 
Eq. (f3"6")l the result is no surprise. At least for small sound 
speeds, gravitational lensing constrains directly the con- 
tribution of the dark energy perturbations to the metric, 
not their impact on the dark matter. 



B. Galaxy power spectrum 

We now explore a second probe of clustering, the 
galaxy power spectrum. Following [29j j we write schemat- 
ically the observed galaxy power spectrum as: 



Pobs\Z, "V) 



n 2 



{Z)H{Z) G\z)b{zf 
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Pshot(z) 



1 



pn 2 yp 0r (k) 

(37) 



where the subscript r refers to the values assumed for 
the reference (or fiducial) cosmological model, i.e. the 
model at which we evaluate the Fisher matrix. Here P s hot 
is the shot noise due to discreteness in the survey, fj, is 
the direction cosine within the survey, P$ r is the present 
spectrum for the fiducial cosmology, G(z) is the linear 
growth factor of the matter perturbations, b(z) is the 
bias factor (assumed scale independent) and Da is the 
angular diameter distance. 

The wavenumber k is also to be transformed between 
the fiducial cosmology and the general one ([29j and 
see also [30j and [3l|, for more details). To avoid non- 
linearity problems (both in the spectrum and in the bias), 
the Fisher matrix is calculated up to a limiting k max (z) 
at z: we choose values from 0.11/i/Mpc for low- z bins to 
0.3/i/Mpc for the highest z-bins. 

Here again the galaxy power spectrum is affected by 
the dark energy perturbations in three different direct 
ways: 

• through the growth factor, see Eq. (j2"Tj) ; 

• through the redshift space distortions, see IIII Dl 

• through the present matter power spectrum P 0r , as 
discussed in IIII CI 

We consider a photometric survey from z = — 2 di- 
vided in equally spaced bins of width Az = 0.2 as our 
benchmark survey; we assume an error on the mea- 
sure of redshift of about Sz/z — 0.01 and an area of 
20000 deg 2 . These features are similar to those of pro- 
posed experiments like JDEM and Euclid, 32[ and [33j j . 
respectively (see also DETF report 20]). Here too we 
also consider extended surveys to z max = 3 and z max = 
4. Moreover, we assume an overall radial distribution 
n(z) = z 2 cxp[— (z/zq) ' ] where zq — z mean /lA12 and 
a z mean = 0.9. The bias factor here is assumed to 
be only redshift dependent; we assume as fiducial bias 
b(z) = s/1 + z. In general, however, dark energy per- 
turbations could lead to a scale dependence of the bias 
factor. Our fiducial model is exactly the same as for the 
WL survey considered in the previous section. 
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Figure 5: Confidence region at 68% for three different value 
of z max = 2, 3, 4, blue dashed, green long-dashed and solid 
counter, respectively. The upper panel shows the confidence 
region when the sound speed is c s = 1; the bottom panel with 
the sound speed c? s = 1CP 5 . The parameter equation of state 
is for both cases u>o = —0.8. 



In Fig. [5] we report the confidence region for w ,c 2 
for two different values of the sound speed and z max . 
For high values of the sound speed (c 2 . — 1) we find, for 
our benchmark survey: <j(w ) = 0.0088, and <j(c 2 )/c 2 — 
22.07. Here again we find that galaxy power spectrum is 
not sensitive to the clustering properties of dark energy 
when the sound speed is of order unity. If we decrease 
the sound speed down to c 2 = 10~ 5 then the errors are 
a(w ) = 0.0091, cr(c 2 )/c 2 = 0.32. 

In Tab. [II] we listed explicitly the errors for Wq and the 
relative errors of c 2 for six different values of the sound 
speed for a benchmark survey up to z max = 2. 

We can ask again the question what part of the power 
spectrum is most sensitive to the dark energy perturba- 
tions, G, (3 or Po(k)l In order to disentangle their relative 
contributions, we plot the corresponding diagonal Fisher 
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Table II: Errors for the equation of state parameter and the 
relative errors for the sound speed for z max — 2. 
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Figure 7: The upper panel shows G — G,f3 — f3,Po — Pq, solid 
blue, red and green lines, respectively. For c 2 = 10 -5 
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Table III: Relative errors for c 2 for P (k) for different z ma x 
without redshift distortion. 



Figure 6: The figures shows the Fisher matrix element G — G, 
P — P, Po — Po, solid blue, red and green lines, respectively. 
For c 2 = 1 



speeds given. For the fiducial surveys we find 

-£■ = 1830 (c?r 



(38) 



matrix terms containing the sound speed as a function 
of redshift. The results are shown in Fig. |6] for c 2 = 1 
and in Fig. [7] for c 2 s = 10 -5 . For large sound speed the 
shape of the dark matter power spectrum is formally the 
most sensitive, but not sufficiently to constrain the per- 
turbations. In the more interesting case where the sound 
speed is lower and the perturbations can be detected, it 
is instead the growth factor that is the most sensitive. 

In both cases the contribution from redshift space dis- 
tortions to the overall constraints is sub-dominant except 
at redshifts below about z — 0.3 and for small c 2 s: cf Fig- 
ures [B] and [7J In the low-redshift region the growth fac- 
tor is less sensitive as it is an integral from z — 0. To 
further clarify the importance of the redshift space dis- 
tortions, we performed our analysis excluding the depen- 
dence of dark energy perturbations on the redshift dis- 
tortion parameter. The result, given in Tab. IIII1 shows 
that for small sound speeds the contribution from /3 is 
non- negligible when we compare it to Table ITT] 



C. The scaling of the errors with c s 

In Fig. [8] the behavior in logarithmic scale of the sound 
speed errors <j{c 2 s )/c 2 s are shown as a function of the sound 
speed itself. We notice that the errors appear to scale as 
a power law for both probes, over the range of sound 



for the weak lensing and 



«C 2 , n,( 



(39) 



for the power spectrum. 



D. Combined Fisher matrices and degeneracies. 

We can also combine the two probes, the galaxy power 
spectrum and weak lensing, to improve the sensitivity to 




Figure 8: The figure shows the relative errors on the sound 
speed for two different experiments and their fits: WL in red 
and P(k) in black, both for our benchmark survey. 
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Table IV: Correlation matrix for galaxy power spectrum ex- 
periments and cl = 10~ 5 . 
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Table V: Correlation matrix for WL experiments and c s 
10" 5 . 



the dark energy perturbations, cf. Tab. IVIII1 As ex- 
pected the errors decrease by about 30% for low values 
of the sound speed and remain unaltered when the sound 
speed is close to unity. 

We find that the sound speed does not suffer of any 
strong degeneracies with the other parameters considered 
in this work. We evaluated the dimcnsionless correlation 
matrix: 



Corr^ = 



Cij 



yj OuG : 



j j 



where Cij is the inverse of the Fisher Matrix Fij 
instance, for c 2 - "' '" 



For 



10" 



we find for all parameters a cor- 
2 of at most 0.3 for galaxy power 



relation factor with c~ 
spectrum experiments and at most 0.23 for WL experi- 
ments; however, for the latter case we find for Q. m< o — c 2 
terms a correlation value of about 0.68. Increasing the 
sound speed the correlation terms decrease for galaxy 
power spectrum case and WL case; for the latter only 



the term ilhh 2 



increases up to 0.78, for cj = 1. 



E. An integrated measure of dark energy clustering 

So far we have been studying the errors on the dark 
energy parameters w and c 2 . A subtly different question 
is whether we are able to detect dark energy perturba- 
tions or not - a question not about parameter constraints 
but about model probabilities [34|]. It is not straight- 
forward to use our constraints e.g. on the sound speed 
for this purpose, as it docs not make sense to consider 
the speed of sound if there are no perturbations at all, 
i.e. a model without perturbations is not nested in the 



quintessence model considered here. This is just another 
way to say that all quintessence models (except those for 
which w — —1, equivalent to a cosmological constant) 
necessarily have perturbations. 

But from the errors on w and c 2 we can forecast the 
constraints on the amount of deviation of the Poisson 
equation from the expectation due to the dark matter 
and other known contributions, given by Q. A model 
without dark energy perturbations would then be char- 
acterized by Q = 1. However, Q is a function which 
evolves and additionally can have a scale dependence. It 
is then useful to define a compact way to express the de- 
viation integrated over z and k. We define therefore the 
average quantity 



w(n m ,w,c 2 s ) 



47T 
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\Q(k,z)- l\dzk 2 dk (41) 



where Az is the observed redshift range and Vu = 
ATrk 2 nax /3 is the momentum volume. Finding a devia- 
tion of W from zero would signal some clustering of dark 
energy. We note that it is straightforward to focus on 
certain scales or epochs by introducing an appropriate 
weight function into the integral above. 

In order to obtain constraints on W we need to 
marginalize the Fisher matrix over all parameters except 
Pi = {fi m , wo, c 2 } and then project it over the new pa- 
rameter set q.j = {W, Wo,c 2 }. Defining Jy = (dpi/dqj) r 
we have 



— (JikFklJli) 



(42) 



The errors on W for the P(k) case are reported in Tab. I VII 
and for WL in Tab. IVIIl Since the derivative of Q with 
respect the sound speed contains a term which goes like 
1/c 2 , the dominating source of error on W comes from 
(40) the errors on c 2 , as long as c 2 <C 1. Therefore a quick 



estimate can be obtained by approximating 
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(43) 



and we find that in the present case the relative errors 
on W are similar to those on c 2 . This supports the naive 
interpretation that we can only detect the perturbations 
if we can measure the sound speed, i.e. if the relative 
error on c 2 is less than unity. 

However, we are now in a position where we can use 
model comparison techniques to decide whether W ^ 
has been detected or not. Using the Savage-Dickey den- 
sity ratio, we find that the relative model probability is 
given by the value of the normalized posterior (marginal- 
ized over all common parameters, i.e. all parameters ex- 
cept W) at W — 0, divided by the value of the normalized 
prior at the same place. To simplify the calculations, we 
choose a uniform prior in all variables except W, and for 
that variable a Gaussian pdf centered at W — with a 
width (variance) E 2 . This then leads to a Bayes factor 
of 



\¥ (w 2 \ 



(44) 



12 



P(k) 


cl 


aw/W 


a^Jcl 


10" 5 


0.15 


0.32 


10- 4 


0.42 


0.66 


10" 3 


1.16 


1.41 


10~ 2 


2.53 


2.79 


lO" 1 


13.75 


14.91 


1 


21.46 


22.07 



Table VI: Relative errors for W and the corresponding c 3 for 
P(k). 
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Table VII: Relative errors for W and the corresponding c s for 
WL. 



Here W is the value on which the posterior for W is cen- 
tered, and a 1 the variance of W. Given the fiducial value 
Wfid and the Fisher matrix error a Wl which characterize 
the likelihood, we find for the posterior the shifted values 
(due to the prior, see e.g. the appendix of Ref. |35j |) 
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However, it is not obvious how we should choose the 
width of the prior, for a general model with unknown 
sound speed. One possibility is to use the most optimistic 
value, for which the model discrimination is maximal [36j. 
To find this value, we maximize Bw with respect to E, 



which leads to E 2 = W| d 



a w . Using this prescription 



(valid only for Wfid > aw) we finally obtain 
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(47) 



This formula depends only on the ratio Wad/aw, and 
if we demand strong evidence in favor of the presence of 
perturbations, In B > 5, we find that we need aw /Wfid < 
0.27, i.e. not quite a 4-ct detection. We remind the reader 
that this corresponds to the choice of prior that max- 
imally favors the detection of perturbations, any other 
choice would need a stronger detection of W to reach the 
same Bayes factor. 

Looking at tables IVI1 1 VIII and IVIIII we can see that in 
our class of quintessence models and for the observations 
that we consider here, we can only hope to strongly favor 
the presence of dark energy perturbations if the sound 
speed is c s < 0.01. 



V. CONCLUSION 

In this paper we extended the idea proposed in [13j . 
We investigated the effect of a generic modification of 
Poisson equation Q(k,a) in the galaxy power spectrum 
and in the weak lensing convergence power spectrum in- 
duced solely by dark energy clustering. In other words, 
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Table VIII: Errors for the equation of state parameter, the 
relative errors for the sound speed and the relative errors for 
the variable W combining WL and galaxy power spectrum 
probes. 



we are using perturbations in the metric and in the dark 
matter to detect the tiny traces of dark energy perturba- 
tions. 

Wc identified several effects: on the power spectrum 
growth, on its shape, and (for the galaxy clustering 
probe) on the redshift distortion. These effect are com- 
pletely characterized by w and c 2 . We performed a Fisher 
matrix analysis of future large-scale surveys to forecast 
the errors on w and c 2 and therefore on Q(k,a) or on 
its average version, W . As perhaps expected, we find 
that dark energy perturbations have a very small effect 
on dark matter clustering unless the sound speed is ex- 
tremely small, c s < 0.01. For c s < 0.01, we find that 
W (that is, the average deviation for a standard Poisson 
equation) could be constrained to within 66% or better. 
Let us remind the reader that in order to boost the ob- 
servable effect, we always assumed w = —0.8: for values 
closer to —1 the sensitivity to c 2 is further reduced; as 
a test we performed the calculation for w 

„2 



-0.9 and 
c; = 10~ 5 and we find c c 2/c 2 = 2.6 and cr c 2/c 2 = 1.09 
for WL and galaxy power spectrum experiments, respec- 
tively. 

Such small sound speeds are not in contrast with 
the fundamental expectation of dark energy being much 
smoother that dark matter: even with c s s=s 0.01, dark 
energy perturbations are more than one order of mag- 
nitude weaker than dark matter ones (at least for the 
class of models investigated here) and safely below non- 
linearity at the present time at all scales. Models of 
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"cold" dark energy with c 2 s — are interesting because 
they can cross the phantom divide [5[ and contribute to 
the cluster masses [HI]. Small c s could be constructed 
for instance with scalar fields with non-standard kinetic 
energy terms. Here we showed that future large scale sur- 
veys have the potential to rule out or confirm this class 
of models. As an example we report here the absolute 
errors for the sound speed (? s = 0: a c 2 = 0.344 • 10 -9 and 
ct c 2 = 0.211 • 10~ 7 for WL and galaxy power spectrum 
experiments respectively. 
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